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Abstract 

Some one-neutron halo nuclei can emit a proton in a /3 decay of the halo 
neutron. The branching ratio towards this rare decay mode is calculated 
within a two-body potential model of the initial core-|-neutron bound state 
and final core+proton scattering states. The decay probability per second 
is evaluated for the ^^Be, ^^C and ^^Ne one-neutron halo nuclei. It is very 
sensitive to the neutron separation energy. 



1. Introduction 

Some neutron-rich halo nuclei can emit a proton. This process is possible 
if the neutron separation energy is very small. Indeed a weakly bound halo 
neutron may (5 decay, producing a proton which can be emitted, possibly 
together with neutrons. Processes where this proton is bound with one or 
two neutrons have been observed in the /3 delayed deuteron and triton decays 
of ^He and ^^Li [H-Isj. Recently we have calculated the branching ratio of an 
even rarer process where the proton remains unbound but is accompanied 
by a free neutron [9']. This decay is uniquely possible for ^^Li, among nuclei 
with known separation energies. The study has been performed in a three- 
body model with a simplified description of the continuum. An even simpler 
process is however possible. 
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A one-neutron halo nucleus can be viewed as a normal nucleus, the core, 
to which a neutron is bound in an orbital with a large radius. The (3 decay of 
the bound halo neutron may occur, releasing the proton, under the condition 
of energy conservation 

< {mn -rup- me)c^ ^ 0.782 MeV, (1) 

where is the neutron separation energy of the decaying nucleus and m„, 
rrip and rrie are the neutron, proton and electron masses, respectively. Among 
one-neutron halo nuclei for which Sn is known with sufficient precision, this 
decay is allowed at least for ^^Be and ^^C, and probably for ^^Na. It should 
be observable if the branching ratio is large enough. This decay mode of ^^Be 
has been considered by Horoi and Zelevinsky but the results do not seem to 



have been published [10[ . Here we study this rare decay mode within a two- 
body potential model. The initial halo nucleus is treated as a core+neutron 
bound state. The final states lie in the core+proton continuum. How rare is 
this decay is the main question raised in the present exploratory study. 



2. Decay probability for f3 delayed proton emission 

The P decay of the halo neutron releases the resulting proton from the 
core. The distribution of decay probability per time unit as a function of the 
energy E < Q oi the relative motion of the two particles is given by 

where Gjb ~ 2.996 x 10~^^ is the dimensionless /3-decay constant and A ~ 
—1.268 is the ratio of the axial-vector to vector coupling constants. The 
Fermi integral f{Q — E) depends on the kinetic energy Q — E available for 
the electron and antineutrino with 

Q = (m„ -rUp- me)c^ - S^- (3) 

The total decay probability per time unit W is obtained by integrating ([2]) 
from zero to Q. The branching ratio can than be derived as 

7^ = ml/2/ln2, (4) 

where ti/2 is the half life of the halo nucleus. 
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In the present model, the halo nucleus is described as a two-body core-l- 
neutron system in its ground state with total angular momentum Jj resulting 
from the coupling of the orbital momentum /j of the relative motion and the 
neutron spin s = 1/2. The spin of the core is assumed to be zero. The parity 
of the initial state is (—1)'*. The radial wave function is denoted as Uil^J^ 
with the normalization \uii^j^{r)\^dr — 1. It is obtained from a potential 
Vi adjusted to reproduce the experimental neutron separation energy S^^. 

The final scattering state of the core and the proton is a distorted wave 
with wave vector k. Because of selection rules, only some partial waves 
with total angular momentum Jf resulting from the coupling of the orbital 
momentum If and the proton spin s are allowed. The radial wave functions 

Ukifjf for a wave number k — ^2fj,E/h'^ where fj, is the core- proton reduced 
mass are obtained with a potential Vf describing the core+proton system. 
They are normalized according to Uki^j^\r)uyi^.j^.{r)dr = S{k — k'). The 
potential Vf is usually poorly known when the core is unstable. 

Within this model, the Fermi reduced decay probability is given by 

and the Gamow-Teller reduced decay probabihty by 

2 



dB{GT) _ 6 ^ r Jf s k \\j |2 



(6) 



with the relative velocity v — hk/ji and the radial integrals 

/■oo 

IiJiJf = / Ukijf{r)uiij^{r)dr. (7) 
Jo 

If the final wave function does not depend on Jf, the Gamow-Teller term 
simplifies as 

^^(GT) dB{F) 

The reduced decay probability can then also be written as 

dW _ f{Q-E) dB{¥) 
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where Wn is the free-neutron (3 decay probabihty per second and /„ is the 
corresponding Fermi integral. 

With respect to a free neutron, the decay probabihty is affected in two 
ways. First, the ratio f{Q — E)/ fn is smaU due to the reduction of phase 
space, since /„ = f{Q + Sn)- It becomes extremely small when E tends to Q. 
The P delayed proton emission is favoured by very small separation energies 
S^. Second, the reduced decay probability appearing in is proportional 
to the square of a radial integral ([7]). Because of the Coulomb repulsion and 
the smallness of the Q value, the scattering waves are small and, when E 
tends to zero, tend to zero as exp(— 7r?7) where r] = ZcC^/hv is the 
Sommerfeld parameter. They become thus smaller with increasing charge Z^. 
of the core. They also become smaller with increasing orbital momentum. 
Hence, at given Q value, we expect the decay probability to be largest for 
the lightest halo nuclei and for the halo neutron in the s wave. 

3. Results and discussion 

Before making explicit calculations, we have to specify the choice of po- 
tentials. The Fermi strength is proportional to the square of an overlap 
integral ([7]) between the initial and final radial wave functions. In order to 
have a realistic overlap, it is useful to have a correct node structure for these 
wave functions. Indeed, the presence of nodes leads to an integrand that 
changes sign one or several times and thus to a reduction of the overlap. 
Spectroscopic factors can also affect the size of the Fermi strength but given 
the limited knowledge on these quantities, we choose to ignore them in the 
present exploratory study. Finally, absorption in the core-|-proton optical 
potential might also play a role. However, the energies of the states after 
decay are lower than, or comparable to, the energy of the Coulomb barrier. 
Absorption should be weak and can safely be neglected. 

Hence, we shall use real potentials Vi and Vj which should be deep enough 
to provide a realistic node structure of the initial and final radial wave func- 
tions. To keep the model simple we only use central Woods-Saxon potentials 
with range roAl^^ where A,, is the mass number of the core. The depth is 
adapted to the separation energy for the core-|-n system. The same form 
factor with an additional point-sphere Coulomb potential is employed for 
the final core+p elastic scattering. Because of the small energies, the phase 
shifts are small and the sensitivity to Vf is weak. Now let us consider explicit 
cases. 
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The best documented case is ^^Be. Its 1/2"'" ground state has a separation 
energy of about 501 keV 12| and its half hfe is 13.8 s 13|. The halo neutron 
is described by an s wave. The parameters of the Woods-Saxon potential 
are taken as tq = 1.2 fm, a = 0.6 fm and Vio = 62.52 MeV ij]. In the 
s wave, this potential possesses one unphysical forbidden state. The same 
parameters are used for the final potential except VfQ. The ^^B nucleus has 



a proton separation energy Sp ~ 11.228 MeV [15|. Its lowest 1/2"*" state is 
located at the excitation energy ~ 6.79 MeV. In the s wave, V/o = 84.1 
MeV is adjusted so that the potential possesses one forbidden state and one 
bound state fitted to the energy — Sp ^ —4.52 MeV with respect to the 
^°Be+p threshold. Bound and scattering states should thus have a reasonable 
node structure. 




The Q value ([3]) is small, 0.281 MeV. The distribution of decay probability 
is displayed in Fig. [T] The most probable energies of the relative motion are 
in the interval 0.1-0.2 MeV. The total decay probability 1.5 x 10""^ s~^ leads 
to a branching ratio 3.0 x 10~^. 
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The ^^C 1/2"'' ground state has a separation energy of 580 ± 90 keV 15 
and a half hfe ti/2 = 46.2 ms 13(|. As a simple picture, we consider a neutron 



in the s wave with one forbidden state and no spectroscopic factor. The 
parameters of the Woods-Saxon potential are tq = 1.25 fm, a = 0.62 fm and 
Vio = 41.42 MeV giving a Q value of 202 keV. For the final ^^C+p system, the 
s wave possesses one forbidden state. We assume a possible 1/2+ bound state 
near = 2.1 MeV With Sp ^ 16.35 MeV 0, we take Vjo = 77.2 



MeV which gives a bound state at —14.2 MeV. The distribution of decay 
probability is displayed in Fig. [H It is much smaller than for ^^Be because 
of the larger charge of the core and the smaller Q value. The total decay 
probability 2.7 x 10^^^ s^^ leads to a branching ratio 1.8 x 10^^^. 

A candidate for delayed proton emission is ^^Ne. Its neutron separation 



energy is poorly known: 0.33 ± 1.07 MeV |15[. Its half life is ti/2 = 3.4 ms 



13|. This nucleus belongs to an island of inversion where its ground state 



should be an intruder state. Its one-neutron removal cross section [17| is too 
large for agreeing with the quantum numbers 0/7/2 of the naive shell model. 
The ground state could be described with a lp3/2 orbital [l8[ although a 
2sl/2 orbital has also been considered j20|. Here we assume a p wave ground 
state at —0.33 MeV giving Q = 0.45 MeV. It can be reproduced with the 



parameters ro = 1.25 fm, a = 0.75 fm and V^o = 48.86 MeV [18|. This 



potential has one forbidden state in the p wave. Little is known about the 
^°Ne+p scattering. One can also expect an intruder 3/2~ state in the vicinity 
of the ground state. Hence we choose V/o = 90.0 MeV which provides a 
forbidden state and a bound state at —16.1 MeV, not far above —Sp ~ —17.7 
MeV. 

The distribution of decay probability is displayed in Fig. [H It is smaller 
than for ^^Be because of the larger charge of the core and the higher orbital 
momentum, but these effects are partly compensated by the larger Q value. 
The most probable energies E lie between 0.25 and 0.35 MeV. The total 
decay probability 3.3 x 10~^° s~^ leads to a branching ratio 1.6 x 10"^^. For 
an s ground state with two forbidden states (Vio = 69.27 MeV), the decay 
probability W ~ 1.6 x 10^^ would be five times larger. 

The separation energy of ^^Ne is quite uncertain. The one-neutron re- 
moval cross section can be interpreted as arising from ~ 0.6 MeV but this 
assumption is weakened by the lack of knowledge of spectroscopic factors 



18| . Hence we display in Fig. [2] the dependence of the decay probability 
on the separation energy S^^. One observes that it varies very strongly. If 
5*11 is around 0.6 MeV, the decay probability is reduced by about six orders 
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Figure 2: Decay probability per second for the (3 delayed np decay of '^^Ne as a function 
of the separation energy S^. 

of magnitude. On the contrary, the decay probabihty can be larger by five 
orders of magnitude if the separation energy is very small. 

Finally, let us note that an estimate of the order of magnitude (in general 
within a factor of two) can be obtained with the simple approximation 



where a = 5 fm, k = y 2/iS'„//i^ and Fi is a regular Coulomb function. Under 
the same conditions as in Fig. [H the asymptotic normalization constant C is 
0.83, 0.96, 0.69 imr^l"^ for ^^Be, ^^C, ^^Ne, respectively. 

4. Conclusion 

As a summary, we have evaluated the order of magnitude of the decay 
probability per second for the (5 delayed proton emission by one-neutron halo 
nuclei. The best candidate for observing such a decay is ^^Be in spite of the 
fact that its separation energy is not very small. The probability of this 
delayed decay is smaller than for the neutron-and-proton delayed decay of 
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^^Li [9] by an order of magnitude. Because of a longer lifetime, the branching 
ratio is larger by two orders of magnitude. The observation of this (3 delayed 
decay mode of ^^Be would thus require high radioactive beam intensities and 
long measurement times. 

The neutron separation energies of the other candidates, ^^C and '^^Ne, 
are less well known and the decay probabilities are thus more uncertain. 
We have shown that the decay probability varies strongly with the neutron 
separation energy. A very small S^i would be advantageous for the study of 
this decay mode. This advantage however decreases when the charge of the 
core increases. The best candidate for observing such a decay would be a 
not too heavy one-neutron halo nucleus with a very small separation energy. 
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